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ABSTRACT 

The  development  of  time  transfer  techniques  using  GPS 
carrier-phase  observations  promises  the  capability  to 
deliver  sub-nanosecond  time  transfer  capabilities. 
NAVSYS’  High-gain  Advanced  GPS  Receiver  (HAGR) 
has  been  designed  to  enhance  the  accuracy  of  the  GPS 
observations  by  using  digital  beam-steering  from  an 
antenna  array.  This  provides  over  10  dB  of  additional 
signal  gain  on  each  satellite  tracked  which  improves  the 
accuracy  of  both  the  code  and  carrier  tracking  loops, 
simplifying  integer  ambiguity  resolution  for  carrier  phase 
time  transfer.  This  same  design  provides  a  highly  phase 


stable  observation  of  the  GPS  carrier,  relative  to  a  local 
reference  oscillator,  that  enables  precise  time  observations 
to  be  made  for  carrier-phase  time  transfer. 

In  this  paper,  the  design  of  NAVSYS’  P(Y)  L1/L2  Code 
HAGR,  developed  for  precise  time  transfer  applications, 
is  described  and  test  results  are  included  showing  the 
accuracy  of  the  code  and  carrier  phase  observations. 

INTRODUCTION 

GPS  carrier-phase  measurements  provide  the  potential  for 
much  improved  precision  in  time  and  frequency  transfer 
[123456].  Because  carrier  phase  observations  are  100  to 
1 000  times  more  accurate  than  code  based  pseudo-range 
measurements,  time  transfer  errors  approaching  100 
picosecond  (ps)  are  expected  using  this  approach. 

The  key  to  reliable  carrier  phase  time  transfer  is  the 
ability  to  resolve  the  carrier  cycle  ambiguity  between  the 
carrier  phase  and  the  pseudo-range  to  the  satellite  and  also 
to  determine  which  carrier  phase  cycle  a  given  receiver  is 
tracking,  related  to  an  external  reference  clock  in  a 
calibrated  fashion. 

This  is  achieved  in  the  P(Y)  High-Gain  Advanced  GPS 
Receiver  (HAGR)  design  through  the  following 
innovations  that  are  described  further  in  this  paper. 

1.  Digital  beam-steering  phased  array.  The  digital 
beam-steering  provides  gain  in  the  direction  of  the 
GPS  satellites,  as  shown  in  Figure  1.  This  increases 
the  measurement  precision,  and  reduces  multipath 
errors  improving  the  ability  to  perform  cycle 
ambiguity  resolution. 

2.  Phase  coherent  Digital  Front-End  (DFE).  The 

phase  coherent  DFE  samples  the  received  GPS  signals 
with  a  clock  precisely  synchronized  and  calibrated  in 
phase  to  an  external  reference.  This  allows 
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unambiguous  phase  observations  to  be  made  coherent 
with  the  external  time  reference  to  be  synchronized. 

3.  Phase  coherent  frequency  synthesis.  The  frequency 
synthesizer  is  designed  to  generate  the  receiver  Local 
Oscillator  (LOs)  for  downconversion  synchronized  to 
the  external  time  reference  both  in  frequency  and  also 
with  an  unambiguous  phase  offset.  This  removes  any 
ambiguous  phase  offset  at  start-up  from  the 
synthesizer  phase-locking  to  the  external  reference. 


inputs  from  an  array  of  antenna  elements.  Digital  signal 
processing  is  then  used  to  create  an  optimal  beam  pattern 
for  each  of  the  satellites  tracked. 

The  receiver  can  be  configured  based  on  the  user’s 
requirements.  In  Figure  2  a  P(Y)  HAGR  configuration  is 
shown  which  is  being  delivered  to  the  US  Naval 
Observatory  for  a  carrier  time  transfer  application.  This 
receiver  is  configured  to  provide  the  following 
capabilities. 


Figure  1  HAGR  Digital  Beamsteering  Array 


•  16-element  digital  beam-steering  on  both  the  LI  and 
L2  frequencies 

•  12  satellite  tracking  channels 

•  C/A  (LI)  and  P(Y)  (LI  and  L2)  observations 

•  Precise  time  synchronization  (code  and  carrier  phase) 
to  an  external  time  reference. 
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Up  to  16  Antenna  6 10  12  Processing 

Elements  Channels 


Figure  2  P(Y)  HAGR  System  Architecture 


P(Y)  HAGR  DESIGN 

The  P(Y)  HAGR  is  based  on  a  modular,  software  and 
firmware  reprogrammable  digital  GPS  receiver  design7. 
The  system  architecture,  shown  in  Figure  2,  uses  an  RF- 
to-Digital  Front  End  (DFE)  to  digitize  each  of  the  antenna 


The  components  of  this  receiver  are  shown  in  Figure  3 
and  are  described  further  below. 
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Figure  3  12-Channel  P(Y)  HAGR  System  Components 


Antenna  Array 

The  antenna  array  consists  of  1 6  dual  frequency  elements 
laid  out  on  a  4x4  antenna  pattern.  The  elements  are 
separated  by  0.5  cycles  at  the  LI  frequency.  As  described 
later,  this  provides  optimal  beam-steering  performance  at 
the  LI  frequency,  while  some  spatial  aliasing  occurs  with 
the  L2  frequency  beam- forming. 


Figure  4  L1/L2  16-element  Antenna  Array 


Digital  Front-End  (DFE)  Cards 

The  DFE  cards  perform  the  RF  to  digital  conversion  on 
each  of  the  antenna  inputs.  Each  card  includes  the 
capability  to  handle  up  to  four  antenna  inputs.  The 


frequency  (LI  or  L2)  of  each  channel  is  selected  by  the 
front-end  filter  and  by  using  either  the  LI  or  L2  LO 
provided  from  the  Local  Time  Generator  (LTG) 
synthesizers.  To  handle  16  RF  inputs  at  LI  and  16  RF 
inputs  at  L2,  the  P(Y)  HAGR  is  configured  with  eight  of 
the  DFE  cards  shown  in  Figure  5. 


Figure  5  Digital  Front  End  (DFE)  Board 


Digital  Beam  Steering  (DBS)  Card 
The  DBS  card  includes  the  digital  logic  to  apply  the 
beam-steering  weights  to  the  LI  and  L2  digitized  RF 
inputs.  Each  card  can  provide  12  digital  beam- steering 
outputs  from  up  to  16  digital  antenna  inputs  from  the  DFE 
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cards.  To  provide  beam-steering  on  12  satellite  channels, 
LI,  and  12  satellite  channels,  L2,  two  DBS  cards  are 
included. 

Correlation  Accelerator  Card  (CAC)  Card 
The  CAC  performs  the  C/A  and  P(Y)  code  generation  and 
correlation  functions  on  the  digital  beam-steered  signal 
outputs  for  each  of  the  satellite  channels.  Each  CAC  card 
is  configured  to  handle  six  satellite  channels.  Each 
channel  provides  the  I/Q  correlator  outputs  (Early,  Prompt 
and  Late)  for  the  C/A  (LI),  P(Y)  LI  and  the  P(Y)  L2 
signals.  P(Y)  code  operation  requires  that  the  CAC  card 
be  keyed  for  PPS-SM  operation  through  a  KYK-13 
interface.  Two  CAC  cards  are  included  to  allow  tracking 
on  up  to  12  satellites  simultaneously. 

Local  Time  Generator  (LTG)  Synthesizers 
In  an  off-the-shelf  RF  phase-locked  synthesizer  design, 
the  local  oscillator  signal  F0  is  phase  locked  to  the  input 
reference  clock  by  feeding  back  to  the  PLL  phase 
comparator  a  divided  version  of  the  output  F0  signal.  In 
commercial  PLL  synthesizer  chips,  this  feedback  division 
is  free  running,  and  as  a  result  the  synthesizer  can  lock  to 
any  number  of  N  ambiguous  initial  phases  (N  is  the 
feedback  division  ratio)  as  shown  in  Figure  6.  Test 
results  showing  this  effect  from  two  LI  synthesizers  is 
shown  in  Figure  7.  Both  synthesizers  are  locked  to  the 
same  10  MHz  reference.  One  of  the  synthesizers  is 
powered  down  and  re-started  periodically.  The  observed 
LO  phase  from  turn-on  to  turn-on  changes  due  to  the 
different  start  times  of  the  synthesizer  divide  chain. 

To  eliminate  this  phase  ambiguity,  NAVSYS  have  come 
up  with  a  proprietary  LI  and  L2  synthesizer  design  that 
forces  a  known  phase  synchronization  of  the  feedback 
divider  chains  of  both  LI  and  L2  LO  channels 
synchronous  to  an  external  event  marker  (typically  a 
1PPS  strobe).  This  removes  the  potential  for  the 
ambiguous  synthesizer  phase  offsets  at  start-up  as  shown 
in  Figure  7. 


Feedback  Divider 


Figure  6  Synthesizer  Frequency  Locked  Divide  Chain 
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Figure  7  Conventional  Synthesizer  turn-on  to  turn-on 
phase  offset 

Host  Computer 

A  COTS  6U  CompactPCI  form  factor  PC  is  used  as  the 
host  computer.  The  system  uses  an  850  MHz  Pentium  III 
CPU  with  up  to  1  Gigabyte  of  DRAM. 

Rack  Assembly 

The  19”  rack  mounted  P(Y)  HAGR  configuration  is 
shown  in  Figure  8.  The  chassis  houses  a  standard  6U 
high  CompactPCI  backplane.  The  bottom  of  the  chassis 
houses  cooling  fans.  The  LTG  synthesizers  are  housed  in 
the  rear  of  the  chassis. 


Figure  8  P(Y)  HAGR  Rack  Mount  Configuration 
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DIGITAL  BEAM-STEERING  OPERATION 


The  digital  beam-steering  is  performed  by  applying 
digital  weights  to  each  of  the  HAGR  antenna  array 
outputs  using  the  DBS  cards.  The  weights  are  applied 
optimized  for  each  satellite  channel,  as  shown  in  Figure  2. 
The  digital  signal  from  each  of  the  HAGR  antenna 
elements  can  be  described  by  the  following  equation. 

Equation  1 

Ns 

yk(t)  =  ^lsi(xk,t)  +  nk(t) 
i= 1 

where  s^x^t)  is  the  ith  GPS  satellite  signal  received  at  the 
kth  antenna  element  and  nk(t)  is  the  noise  introduced  by 
the  kth  DFE  at  the  kth  antenna  element 


The  GPS  satellite  signal  at  each  antenna  element  (x^)  can 
be  calculated  from  the  following  equation. 


Equation  2 

2n  t 

Si(xk,t)  =  s,  (0,  t)  exp  {-/—!,  Xk)  =  si(0,  t)esik 

A 


where  Si(0,t)  is  the  satellite  signal  at  the  array  center  and 
li  is  the  line-of-sight  to  that  satellite 
esik  are  the  elements  of  a  vector  of  phase  angle 
offsets  for  satellite  i  to  each  element  k 


The  combined  digital  array  signal,  z(t),  is  generated  from 
summing  the  weighted  individual  filtered  DFE  signals. 
This  can  be  expressed  as  the  following  equation. 


Equation  3 

'  Ns 

^Si(t)esi  +n(t) 

.  i= 1 

With  beam- steering,  the  optimal  weights  are  selected  to 
maximize  the  signal/noise  ratio  to  the  particular  satellite 
being  tracked.  These  are  computed  from  the  satellite 
phase  angle  offsets  as  shown  in  the  following  equation. 


t  t 

z(t)  =  w  y(t)  =  w 


Equation  4 


exp|-'T',r ~]) 

expH-yl  ]xM) 


In  Figure  9  and  Figure  10,  typical  LI  and  L2  antenna 
patterns  are  shown  for  four  satellites  in  view  using  the  1 6- 
element  beam-forming  antenna  array.  As  the  satellite 
move  across  the  sky,  the  beams  follow  the  satellites 
dynamically  adjusting  the  gain  on  each  of  the  channels  to 
optimize  the  signal/noise  ratio  based  on  that  satellite 
direction. 


Figure  10  L2  Beam-Steering  Pattern 

CARRIER  PHASE  TIME  TRANSFER 
PERFORMANCE 

The  GPS  LI  and  L2  pseudo-range  and  carrier-phase 
observations  are  described  by  the  following  equations. 

Equation  5 

PRa(m)  =  Rf  +bu  +  SbuX  -bn  +It  +An  +rM\i  +nPR\ 

f 2 

PRn(m)  =  Rt  +K  +Sbu 2  -bn  -5ba  +It  -^-  +  ATi  +rM2i  +nPR2 

J  2 

CPHn  (m)  ~NXAX+  nCPm  - 

(Ri  +K  +SbcX  -bn  -Sbicl  -Ii  +ATi  +Ax6Mli) 
CPHi2(m)  =  N2 

^2  +  nCPH2 

f2 

(Ri  +bu  +Sbc 2  -biX  - Sbic2  -/;  — ^-  +  Ar-  +A20M2i) 

J  2 

The  following  errors  affect  the  pseudo-range  and  carrier 
phase  observations. 
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1.  Ionosphere  errors-  these  are  different  on  the  LI  and 
L2  PR  and  CPH  observations  (I  and  /  /,2  / /22  ) 

2.  Troposphere  errors  -  these  are  the  same  on  all  of  the 
observations  ( ATi ) 

3.  Receiver  Measurement  Noise  -  these  are  different  on 
each  of  the  observations 

(  nPR\  I  nPR2  I  wCP/fl  »  nCPH2  ) 

4.  Multipath  Noise  -  these  are  different  on  each  of  the 
observations  ( rMli ,  rM2i ,  Ax 0Mli ,  h20M2i  ) 

5.  Satellite  and  Station  Position  error  -  these  affect  the 
ability  to  correct  for  the  Range  to  the  satellite  (Ri) 

6.  Station  clock  and  phase  offsets  on  the  LI  and  L2 
signals  -  these  will  cancel  in  the  time  transfer 
application  (  bn ,  Sbi2 ,  Sbicl ,  Sbic2 ) 

7.  Receiver  clock  calibration  errors  between  the 
reference  clock  offset  (bu)  and  the  LI  and  L2  PR  and 
CPH  observations  ( Sbul ,  Sbu2  ,  Sbcl ,  Sbc2  ) 

Test  results  are  described  in  the  following  sections  that 
characterize  the  clock  calibration  stability,  and  the 
receiver  measurement  and  multipath  noise  from  the  P(Y) 
HAGR  digital  beam- steering  GPS  receiver. 

CARRIER  PHASE  OBSERVATION  ACCURACY 
AND  PHASE  CALIBRATION  STABILITY 

The  accuracy  of  the  carrier  phase  observations  and  the 
clock  calibration  stability  was  observed  by  connecting 
two  receivers  to  the  same  GPS  antenna  and  the  same 
external  time  reference.  The  carrier  phase  observations 
between  the  receivers  were  then  differenced  to  observe 
the  receiver’s  carrier  phase  measurement  error  and  any 
changes  in  the  calibration  parameters  between  the  two 
receivers.  The  results  of  these  tests  for  the  different 
satellites  tracked  are  shown  in  Figure  11  and  Figure  12 
and  summarized  in  Table  1.  This  testing  shows  that  the 
carrier  phase  observations  provided  by  the  HAGR  can 
support  picosecond  time  transfer  resolution. 


Table  1  Summary  of  Carrier  Phase  Time  Transfer 
Observations 
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Figure  11  Carrier  Phase  Difference  between  Two 
Units  (Test  1)  (cycles) 
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Figure  12  Carrier  Phase  Difference  between  Two 
Units  (Test  2)  (cycles) 

PSEUDO-RANGE  MEASUREMENT  NOISE  AND 
MULTIPATH  ERRORS 

The  accuracy  of  the  HAGR  pseudo-range  observations  is 
a  function  of  the  received  signal  strength.  A  data  set  was 
collected  to  observe  the  signal-to-noise  ratio  on  the  C/A 
and  P(Y)  code  HAGR  data  over  a  period  of  12  hours. 
From  this  data  (Figure  17  and  Figure  18)  it  can  be  seen 
that  the  beam-steering  increases  the  GPS  signal  strength 
to  a  value  of  56  dB-Hz  on  the  C/A  code.  As  expected  the 
P(Y)  code  observed  signal  strength  is  3  dB  lower.  The 
predicted  pseudo-range  noise  expected  at  these  signal 
strength  levels  is  shown  in  Figure  15.  The  test  data  was 
analyzed  to  observe  the  pseudo-range  noise  and  compare 
it  against  these  predicted  accuracies. 

From  Equation  5,  the  LI  pseudo-range  +  carrier  phase 
sum  cancels  out  the  common  errors  and  the  range  to  the 
satellite  and  observes  the  pseudo-range  and  multipath 
errors  as  well  as  the  change  in  the  ionospheric  offset. 


6 


PRn  +  CPH  a  ( m )  —  Sbul  +  2/.  +  rM1.  +  +  NxXx  + 

-(^cl  ~~  ®°>ic\  +A  dM\i) 

=  C  +  2/.  +tmu  +nPRX  +(nCPm  ~  @mu  ) 
~  C  +  21.  +  rM1.  + 


The  PR+CPH  is  plotted  in  Figure  16  for  each  of  the 
receiver  data  sets.  The  short  term  (<100  sec)  white 
receiver  noise  was  removed  by  passing  the  PR+CPH 
observation  through  a  linear  filter.  The  drift  caused  by 
the  ionosphere  on  each  observation  was  removed  using  a 
polynomial  estimator.  The  remaining  cyclic  error  is  an 
estimate  of  the  multipath  pseudo-range  errors.  The  peak- 
to-peak  cyclic  PR  variation  for  each  of  the  receiver  data 
sets  was  calculated  used  to  estimate  the  errors  observed 
for  each  satellite  from  the  pseudo-range  multipath8. 
These  errors  are  listed  in  Table  2  for  each  of  the  receiver 
data  sets. 

The  RMS  white  noise  on  the  pseudo-range  observations 
was  computed  by  differencing  the  PR+CPH 
measurement.  This  is  shown  in  Figure  17  and  Figure  18 
for  all  of  the  satellites  tracked  for  the  C/A  and  P(Y)  code 
observations.  The  observed  PR  noise  shows  good 
correspondence  with  the  predicted  values  shown  in  Figure 
15.  For  C/NO  values  above  52  dB-Hz,  the  P(Y)  code 
HAGR  provided  pseudo-range  accuracies  of  5  cm  (1- 
sigma)  while  for  C/NO  values  above  55  dB-Hz  the  C/A 
code  observations  were  accurate  to  15  cm.  These  values 
are  for  1-Hz  observations  without  any  carrier  smoothing 
applied.  The  mean  observed  RMS  accuracies  are 
summarized  below  in  Table  2  with  the  average  peak 
multipath  PR  errors  observed. 


C/A  HAGR 


Figure  13  C/A  HAGR  Signal-to-Noise  (dB-Hz) 


Time  (hrs) 

Figure  14  P(Y)  HAGR  Signal-to-Noise  (dB-Hz) 


Figure  15  C/A  and  P(Y)  HAGR  RMS  PR  error  versus 
C/NO 
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Figure  16  PR+CPH  (m)  -  SV  25 
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CONCLUSION 


Time  (hrs) 


Figure  17  HAGR  C/A  Code  Pseudo-Range  Noise  (m) 
(1-Hz  DLL  -  no  carrier  smothing ) 


P(Y)  RMS  PR  Noise  (m) 
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Figure  18  HAGR  P(Y)  Code  Pseudo-Range  Noise  (m) 
(1-Hz  DLL  -  no  carrier  smothing ) 


Table  2  Mean  PR  Noise  and  M-path  Peak  Errors  (m) 


SVID 

C/A 
HAGR 
RMS  PR 

C/A 
Mean 
Mpath  PR 

P(Y) 
HAGR 
RMS  PR 

POO 

Mean 
Mpath  PR 

1 

0.239 

0.259 

0.054 

0.202 

3 

0.284 

0.494 

0.056 

0.337 

8 

0.200 

0.278 

0.045 

0.202 

11 

0.278 

0.535 

0.059 

0.287 

13 

0.252 

0.321 

0.059 

0.260 

14 

0.214 

0.359 

0.049 

0.350 

20 

0.222 

0.267 

0.050 

0.164 

21 

0.252 

0.261 

0.058 

0.133 

22 

0.248 

0.318 

0.047 

0.217 

25 

0.202 

0.362 

0.044 

0.265 

27 

0.183 

0.270 

0.044 

0.178 

28 

0.236 

0.366 

0.055 

0.272 

29 

0.225 

0.312 

0.050 

0.217 

30 

0.477 

0.791 

0.089 

0.624 

31 

0.325 

0.266 

0.055 

0.135 

The  testing  performed  to  date  on  the  P(Y)  HAGR  digital 
beam-steering  GPS  receiver  has  demonstrated  the 
following  capabilities. 

•  Measure  the  GPS  carrier  phase  noise  to  an 
accuracy  of  6  mm  1 -sigma  (20  pico-seconds). 

•  Measure  the  GPS  P(Y)  pseudo-range 
observations  to  an  accuracy  of  5  cm  1 -sigma 
(166  pico-seconds)  when  the  signal/noise  ratio  is 
above  52  dB-Hz  (includes  +10  dB  gain  from 
beam-steering) 

•  Reduce  the  effect  of  code  multipath  through 
directional  beam-steering,  which  will  improve 
the  ability  to  perform  carrier  cycle  ambiguity 
resolution 

The  combination  of  these  capabilities  is  expected  to 
improve  the  robustness  and  precision  of  GPS  carrier  time 
transfer  by  improving  the  GPS  measurement  quality  and 
the  ability  to  perform  carrier  cycle  ambiguity  resolution. 
The  HAGR  is  also  designed  to  provide  precise  time  and 
phase  synchronization  to  an  external  time  reference. 
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